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Abstract 

An event sample enriched in semileptonic decays of b hadrons is selected using an 
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^ ■ inclusive lepton selection from approximately 3.0 million hadronic Z'^ decays collected 

' ' with the OPAL detector. This sample is used to investigate B meson oscillations by 

Qj I reconstructing a proper decay time for the parent of each lepton, using a jet charge method 

• to estimate the production flavour of this parent, and using the lepton charge to tag the 

decay flavour. We measure the mass difference between the two B[| mass eigenstates 



Amd = 0.444 ± 0.029 tomr ps"^ • 

For the Eg system, we find Arris > 3.1 ps~^ at the 95% confidence level. This limit varies 
only a little if alternative limit setting approaches are adopted. Regions at higher Am-s 
values are also excluded with some methods for setting the limit. 

By studying the charge symmetry of the B[} mixing structure, we are able to constrain 
possible CP and CPT violating effects. We measure the CP violation parameter 

Re €B = -0.006 ± 0.010 ± 0.006 

and the indirect CPT violating parameter 

Im^B = -0.020 ± 0.016 ± 0.006 . 

If we invoke CPT symmetry, then we obtain 

ReeB = 0.002 ± 0.007 ± 0.003 . 
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1 Introduction 



The phenomenon of B — B mixing is now well established. Particle-antiparticle oscillations 
arise when the weak eigenstates, |B°)Q and |B°), differ from the mass eigenstates, |Bi) and 
IB2), which can be described as follows: 

(l + e + 5)|B0) + (l-e-5)|BO) 



|Bi) 



2(1 + |e + 5| 



2^ 



IB2) = vl + ^-'^)|B°)-(l-e + '^)|B°) 
/2(l + |e-5|2) 



where e and 5 are complex, and parametrise indirect CP and CPT violation Note that e 
allows for CP and T violation while respecting CPT symmetry, and 6 allows for CP and CPT 
violation, but respects T symmetry. This formalism applies to both the Bj and B^ systems, 
with separate values for e and 6 in each system. The mass difference, Arriq, between the two 
mass eigenstates for the Bq-Bq system determines the frequency of oscillation, where Bq stands 
for either B° or B^. In the Standard Model, transitions between B^ and 6° mesons arise 
dominantly via box diagrams involving virtual top quarks. Predictions for the mass differences 
can be made, and depend on the top mass and the CKM elements Vtq. For the B° system, 

Amd oc floB^omlF{m^,)\VtiV,l\' (2) 

d cl 

where the first two factors are the meson decay constant and QCD bag model vacuum insertion 
parameter, respectively. These factors are obtained from lattice QCD calculations and QCD 
sum rules, but there is an uncertainty of the order of 50% on this product ||2|. The next two 
factors are the top quark mass, rrit, and a known quadratic function of rrit. The last factor is 
a product of CKM matrix elements. At this point, an accurate measurement of Amd alone 
would not lead to an accurate result for Vtd because of the uncertainty in the decay and bag 
constants. However, if both Am-d and Arris could be measured then information on the CKM 
matrix could be extracted via the ratio 

Arris ^ rriB2_ . l^sP _ /bQ-^bQ 
Amd rriBO iKdP /ro^bo ' 

where m^o and mgo are the B^ and B° masses, as the ratio of decay constants for B[J and B^ 
mesons is much better known than the absolute values [0, 

CP violation has so far been observed only in the K° system, but it is also expected to 
occur in the B° system. Predictions for Reee are of the order of 10~^ in the Standard Model 
for the B° system |^ , and up to an order of magnitude larger in the superweak model . The 
analogous parameter for the Bg system. Re eeo , is expected to be smaller by at least an order 
of magnitude, because Im A ■ Im V^d, where A comes from the Wolfenstein parametrisation 
of the CKM matrix [Q. CPT violation may occur in certain string description models of 
fundamental particles Although CPT conservation has been tested with high precision in 
the K*^ system 0, CPT violating effects could be larger in the B system ^. 

^ We use B° to refer to cither of the neutral B mesons, B[J and Bg. 
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The study presented in this paper uses inclusive lepton events, where the identified lepton 
serves both to select Z — bb events and to determine the flavour of the parent b hadron at 
decay. The production flavour of this b hadron is determined using a jet charge technique. A 
vertex flnding algorithm is used to estimate the decay length of the b hadron, and a separate 
algorithm is used to estimate the energy of the hadron, allowing the decay proper time to be 
calculated. The study of Arrid and Arris is performed neglecting possible CP or CPT violation. 
Results on Reee and Im^B for the system are then obtained, assuming a negligible CP 
violating contribution from decays. 



2 The OPAL Detector 



The OPAL detector has been described elsewhere 0, |ri|]. Tracking of charged particles is 
performed by a central detector, consisting of a silicon microvertex detector, a vertex chamber, 
a jet chamber and z-chambers.0 The central detector is positioned inside a solenoid, which 
provides a uniform magnetic fleld of 0.435 T. The silicon microvertex detector consists of two 
layers of silicon strip detectors; the inner layer covers a polar angle range of | cos^| < 0.83 and 
the outer layer covers |cos^| < 0.77. This detector provided both 0- and ^-coordinates for 
data taken in 1993 and 1994, but 0-coordinates only for 1991 and 1992. Only 0-coordinate 
information was used in this analysis. The vertex chamber is a precision drift chamber which 
covers the range | cos^| < 0.95. The jet chamber is a large-volume drift chamber, 4.0 m long 
and 3.7 m in diameter, providing both tracking and ionisation energy loss (dE/dx) information. 
The 2;-chambers measure the 2;-coordinate of tracks as they leave the jet chamber in the range 
|cos6'| < 0.72. The coil is surrounded by a time-of-flight counter array and a barrel lead- 
glass electromagnetic calorimeter with a presampler. Including the endcap electromagnetic 
calorimeters, the lead-glass blocks cover the range |cos0| < 0.98. The magnet return yoke 
is instrumented with streamer tubes and serves as a hadron calorimeter. Outside the hadron 
calorimeter are muon chambers, which cover 93% of the full solid angle. 



3 Event simulation 



Monte Carlo events are used to predict the relative abundances and decay time distributions for 
lepton candidates from various physics processes. The JETSET 7.4 Monte Carlo program [|12 

was used to generate Z° 



with parameters tuned to OPAL data |p!3| 
subsequently processed by the detector simulation program 



qq events which were 
In this version, branching 

The 



15 



fractions of heavy hadron decays were revised better to reflect measured results 
fragmentation of b and c quarks was parametrised using the fragmentation function of Peterson 
et al. |T6[, with (xe) for weakly-decaying b and c hadrons given by the central values in Table |1|. 



Standard Model values of the partial widths of the Z° into qq were used The mixture 
of c hadrons produced both in Z° cc events and in b-hadron decays was as prescribed in 



^ The coordinate system is defined with positive z along the e beam direction, 9 and 
azimuthal angles. The origin is taken to be the centre of the detector. 



being the polar and 
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Quantity 


Value 


{xE)h 


0.701 ±0.008 [17 


{xe)c 


0.51 ±0.02 ig 


B(h i) 


(10.5 ±0.6 ±0.5)% in 


B{h^c^ i) 


(7.7 ±0.4 ±0.7)% ig 


B{h^c^ £) 


(1.3 ±0.5)% |T|,|li 



Table 1: The parameters used for the Monte Carlo simulation. The uncertainties quoted 
represent the measurement errors, as found in the quoted references. 



reference |]T8[ . The semileptonic branching ratios of charm hadrons and associated uncertainties 



were also those of reference |T^. The central values in Table |I| were taken for the inclusive 
branching ratios for h i, h ^ c —>■ i and b ^ c ^ ^. The models used in describing the 
semileptonic decays of b and c hadrons were those used in determining the central values in 
reference |jl8[. Both B** and D** production were included in the simulation, with production 
rates of 36% of inclusive B and D meson production. Although the measured values quoted in 
Table ^are not in all cases the most recent, the differences are not significant for this analysis. 

Track impact parameter resolutions were degraded to bring the distributions into agreement 
with data, as described previously 



4 Event selection 

The analysis was performed on data collected in the vicinity of the Z° peak from 1991 to 1994. 
Hadronic Z° decays were selected using criteria described in a previous publication . For this 
analysis, events were rejected if the silicon microvertex detector was not operational, resulting 
in a sample of about 3.0 million events. Tracks and electromagnetic clusters not associated to 
tracks were grouped into jets using a cone jet algorithm |2^. The size of the cone was chosen 
so as to include nearly all the decay products of a b hadron into one jet. The jets also include 
some particles produced in the fragmentation process, which originate directly from the e^e~ 
collision point. 



Electrons and muons with p > 2 GeV/c were identified as described in reference A 
sample enriched in semileptonic decays of b hadrons was selected by requiring akin > 0.7, where 
akin is the output of a neural network based on kinematic variables EM. The inputs to this 



network were p, pt and a measure of the energy around the lepton, where pt is the transverse 
momentum relative to the jet axis (calculated including the lepton). The only difference com- 
pared to the previous paper is that the requirement | cos^| < 0.9 was removed. Of the final 



event sample, just over 3% of the leptons have | cos^| > 0.9, ranging up to | cos 6*1 = 0.96. 

For each lepton, an estimated proper decay time was reconstructed for the supposed parent 
b-hadron as described previously |Q. To summarise, a secondary vertex is reconstructed in 



the X — y plane using an algorithm to combine tracks with the lepton track. The algorithm 
was able to form secondary vertices for 70% of the lepton candidates; the remaining 30% 
were discarded. A decay length is reconstructed by fitting to the reconstructed primary and 
secondary vertex positions, using the jet 3jX1S clS db direction constraint. The 3-dimensional 
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decay length is calculated by dividing the 2-dimensional decay length, L^y, by the sin6' of 
the jet axis. The energy of the b hadron is estimated by first reconstructing the energy of 
the jet that includes the lepton, using the Z° mass to constrain the event kinematics, and 
then subtracting an estimated contribution from fragmentation particles. The fragmentation 
particles were separated from the b-hadron decay products using momentum, angle and vertex 
information. The proper time, t, is formed from the (3-dimensional) decay length L and boost |: 

t = -^= , L. (4) 

01 J El - ml 



In this analysis, use is also made of the estimated uncertainty, cxt, on the proper time, calculated 
from the separately estimated uncertainties on the decay length, ctl, and the boost factor, ap^: 




(5) 



where correlations between the uncertainties on L and /?7 are neglected. This is unimportant 
because the shape of the t distribution is parametrised from Monte Carlo, as discussed in the 
next section. The use of at-, estimated event by event, improves the sensitivity of the analysis 
since only events with small at can discriminate between large oscillation frequencies. For each 
event, an initial estimate of ai is calculated from the error matrices of the tracks assigned to 
the secondary vertex. Monte Carlo studies showed that the decay length resolution worsens as 
the total momentum of tracks assigned to the vertex (excluding the lepton) decreases. These 
studies also showed that the resolution worsens with increasing L.j.y. Both effects result from the 
inclusion in the secondary vertex of tracks that do not originate from the decay of the b hadron 
or its daughters, i.e. fragmentation products. Neither of these effects is fully described by the 
initial estimate of a^. Corrections were therefore applied to a^ based on parametrisations of the 
Monte Carlo predictions. Monte Carlo studies of the b-hadron energy estimate showed that the 
resolution depends strongly on L^y and the fragmentation energy in the jet. The uncertainty 
CT/j-y was estimated by parametrising the dependence on these two quantities predicted by the 
Monte Carlo. The performance of the proper time reconstruction is shown in Figure |^ for 
Monte Carlo events, and the distribution of at is shown in Figure ^ 

For each event, the b-hadron flavour at production was tagged using the jet charge technique 



described in a previous paper [25]. The jet charge was defined as 



Q^et = Y.^^{^] , (6) 

where g^, p\ are the charge and longitudinal momentum component (along the jet direction) 
of track z, -Ebeam is the beam energy, and the sum is over all tracks in the jet. The jet charge 
was formed for two jets: the jet containing the lepton and the highest energy jet that did not 
contain the lepton, with n chosen to be and 1 respectively. These two values of jet charge 
were combined into a single value, 

Q2,et = - 10 X Q]:t\opp) , (7) 

where £ denotes the lepton jet, and 'opp' denotes the other jet. The factor of 10 was chosen to 
give a charge separation between produced b and b hadrons that is close to optimal 

^ We use the notation Ti — c — 1. 
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OPAL Monte Carlo 

t' < 1 ps 2 ps < t' < 4 ps 




-10 -5 5 10 -10 -5 5 10 




(t-t')/ps (t-t')/ps 

Figure 1: Distributions of t, (t — t')/o't and t ^ t' for scmilcptonic decays of b hadrons in two 
ranges of true decay proper times: t' < 1 ps and 2 ps < t' < 4 ps. The points represent the 
Monte Carlo simulation, and the curves represent the parametrisation. 
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- — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — - 




a^/ps 

Figure 2: Distribution of at- The points represent the data, and the hne the Monte Carlo 
prediction. 



If more than one lepton was selected in a given event, only the one with the highest value of 
( j^)^ +pf was retained, where p is in GcV/c. A total of 94 843 events remain after all selection 
criteria, where about 75% of the lepton candidates are expected to originate from semileptonic 
decays of b hadrons. 



5 Likelihood function 

The results obtained in this paper are extracted using a maximum likelihood fitting procedure, 
where the overall likelihood is the product of the likelihoods calculated for each event. In order 
to construct the likelihood function, the contributions to the selected lepton candidates were 
split into 11 sources: 

1. lepton candidates from light quark events (Z — > uu, dd and ss); 

2. leptons from semileptonic decays of c hadrons in Z — > cc events; 

3. lepton candidates in Z ^ cc events not from semileptonic decays of c hadrons; 

4. lepton candidates in Z — > bb events that do not originate from a b hadron or its daughters; 

5. lepton candidates in Z — bb events that do not come from semileptonic decays of b 
hadrons or c hadrons, but do originate from a b hadron or its daughters; 

6. leptons from b — > c — > £ decays originating from decays; 
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7. 


leptons 


from b ^ c — s> £ decays originating from 3° decays; 


8. 


leptons 


from b ^ c — s> £ decays originating from B"*' or b-baryon decays; 


9. 


leptons 


from semileptonic B° decays; 


10. 


leptons 


from semileptonic decays; 


11. 


leptons 


from semileptonic B+ or b-baryon decays. 



Note that the notation h —>■ c ^ i used here excludes decays of the type h —>■ c ^ i, where the 
b refers to the quark that decays, after any mixing has occurred. The last three sources are 
taken also to include decays of the type b ^ c ^ £ as well as b ^ r — £. Note also that decays 
of the type h —>■ J ftp ^ £ are classed either as b ^ c ^ £ or semileptonic b hadron decays 
depending on the charge correlation of the lepton and the b quark. The predicted fraction of 
the overall sample contributed by each source is indicated in Table 0. The statistical errors on 
these numbers are not important. 



Source 


1 


2 


3 


4 


5 


6 


7 


8 9 


10 


11 


Fraction 


0.063 


0.089 


0.011 


0.005 


0.013 


0.025 


0.005 


0.025 0.299 


0.089 


0.377 



Table 2: The predicted fraction of the overall selected sample due to each source. 



The joint distribution of t and at is of the form 

%,{t,at) = J S,{at\t')nMc^t,t')V,j{t%{t')dt' (8) 

where t' is the true proper decay time for a given source, V is the physics function, S is an 
efficiency function, TZ is the resolution function for given values of t' and at, and S gives the 
distribution of at for a given t'. The subscript i represents the source number (1 to 11). When 
CP violation is neglected, the subscript j is 2 for leptons from decays of B mesons that decayed 
as the antiparticle of the meson produced (mixed), or 1 otherwise. Note that the physics 
function describes the evolution of each source as a function of the true proper decay time. 
For sources 1, 3 and 4, V is taken to be a delta function at t' = 0. Any true lifetime content 
in these sources is absorbed by the resolution function. For source 2, t' represents the proper 
decay time of the c hadron, but for all the other sources t' represents the proper decay time 
of the b hadron. Sources 6,7,9 and 10 are divided into two components: unmixed (j = 1) and 
mixed (j = 2). For these components: 

_ exp(-^) 1 + cos Am/ _ exp(-^) 1 - cos Am/ 

n ' 2 ' r, ■ 2 ' 

where Amq = Amj (Amg), and Xj = rgo (tbo) for sources 6 and 9 (7 and 10). When consid- 
ering possible CP and CPT violation, sources 6 and 9 (B^) were divided into four components 
according to the flavour of the B at production and decay M : 



^0 uo ^ ^M-tJ l + cosAmdt'-4Im(5BsinAmdt' 
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exp(-^) l-cosArridt' l-4ReeB 

Ti 2 kB 

exp ( - ^ ) 1 + cos Am^t' + 4Im 5b sin Am^t' 

Ti 2/cb 
exp(-^) 1-cosAmdt' l + 4ReeB 
2 

where is a constant which ensures the normahsation J{Vii + Vi2)dt' = 1, and similarly 
for fcg and Vi^ + V^. Note that CP and CPT violation in the Bg system is neglected in this 
analysis. Sources 8 and 11 were described by a double exponential, with the lifetime parameters 
corresponding to the B"*" and b-baryon lifetimes. Source 2 was also parametrised by a double 
exponential to describe the different individual lifetimes of the c hadrons. This approximation 
gave a good description of this background in the Monte Carlo. 

The function £ describes the efficiency for reconstructing a proper decay time as a function 
of t' . A single function was common for sources 9,10,11, a second function was common for 
sources 6,7,8, and individual functions were used for sources 2 and 5. For sources 9,10,11, the 
efficiency is about 5% smaller at t' = ps than at t' = 1.5 ps, above which point it is flat. 

The functions S and TZ were parametrised using Monte Carlo events. For this purpose, 
one set of functions was common to sources 5,9,10,11, and another set was common to sources 
6,7,8. The effects of leptons from secondary decays (b— s>c^£, h ^ c ^ ^, h ^ t tj are 
taken into account in the resolution functions, rather than the physics functions. Distributions 
of t, t — t' and {t — t')/at for semileptonic decays of b hadrons are shown in Figure |I] for two 
ranges of t' . These are compared to the resolution function IZ integrated over the relevant 
proper time interval, and also over at- The distribution of at is shown for data and Monte 
Carlo in Figure ^, equivalent to S integrated over t' and averaged over all sources. The Monte 
Carlo describes the data well, and it is clear that some events are measured much better than 
others. 

The normalisation of the joint distribution function T for each source is determined by 

I J2Vi,it')£,{t')dt' = 1 , lniit\at,t')dt = l , j Si{at\t')dat = I . (11) 

For sources involving mixing, the normalisation condition for V applies to the sum over mixed 
and unmixed functions. 

To describe the charge distributions for each source, we flrst deflne the fraction of each 
source, ^jj, for which the lepton has the same charge sign as the quark, at its production, that 
produced the lepton jet. For example ^ is 1 for semileptonic decays of b hadrons, but for 
b — s> c ^ decays when no mixing has occurred in either case. The full set of values for ^ 
is given in Table ^, where the values for sources 1,3,4 and 5 are taken from the Monte Carlo 
simulation. Variation of these values was found to lead to negligible changes in the results of 
this paper. The sign of the jet charge, Q2jet, is an estimate of the charge sign of the quark, 
before any mixing, that produced the lepton jet. The quantity ?7ij ( | Q2jet | ) is deflned as the 
probability that these charge signs disagree, where i represents the source number and j = 1 
or 2 (also 3 or 4 when considering CP(T) violation as in equation 10) to represent no mixing 
or mixing, respectively, in the lepton jet. Monte Carlo events were used to parametrise r/ as a 



B^^BO : = 
BO^B" : P,3 = 
BO->B° : = 
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function of \Q2jet\- For sources from Z ^ bb events, i.e. sources 4 to 11, the effect of mixing in 
the jet that does not contain the lepton was taken into account by modifying rj as a function of 
X, the time-averaged mixing parameter, averaged over all b hadron species weighted by their 
production fractions in Z" decays: 

T]^ = r,f=' . (1 - X) + (1 - r/,f°) ■ X , where X = x ■ w{\Q2,et\) , (12) 

and the weighting function w{\Q2jet\) was parametrised from Monte Carlo data. Note that x 
is computed in the fit from Amd, Amg, the lifetimes and the fractions of B[j and Bg mesons 
produced. To describe the charge distribution for each source, we define 

Q . = [ ^iyO-- Vij) + (1 - ^ij) ■ Vij if QeQ2jet > .^gs 

\ {1 - ^ij) ■ (l - Vij) + iijVij otherwise 

where is the charge of the lepton candidate. The function Qij represents the probabilities 
for finding QiQ2jet positive or negative for a given \Q2jet\- A comparison of the distribution of 
QeQ2jet between data and Monte Carlo is shown in Figure 0, where the Monte Carlo distribution 
assumed Amd = 0.45 ps~^ and Arris = 15 ps~^. A reasonable agreement is observed. The 
distributions of QeQ2jet for semileptonic decays of neutral B mesons that have and have not 
mixed are compared in Figure^. For neutral unmixed B mesons, 7] is about 0.3 for \Q2jet\ = 2. 
When considering CP(T) violation, two out of the four physics functions defined in equation 10 
are used for a particular source, according to the lepton charge. For example, for a negatively 
charged lepton, the functions corresponding to j = 1 and j = 4 would be used for source 8. 



Source 


1 2 3 4 5 6789 10 11 


^ unmixed 
mixed 


0.69 1 0.62 0.67 0.67 1 1 1 
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Table 3: Values of ^ for each source. The second row represents the case of leptons from B 
mesons that have mixed. 



The final ingredient to the likelihood is the fraction of each source, which is estimated event 
by event based on the value of the neural net output, ctkin, with a small correction for the 
absolute value of the jet charge, \Q2jet\, giving 

The function Ai is the fraction of events due to source i for a given value of akin, and Bi describes 
the distribution |Q2jet| for source i. Note that different sources give rise to slightly different 
distributions of \Q2jet\, in particular light quark events give rise to larger values of |(52jet| than 
do bb events on average. The correction, using the function Bi, takes this into account. Sources 
9, 10 and 11 give rise to very similar distributions of \Q2jet\, and are described by a common 
function B. The function ^i(akin) is parametrised from Monte Carlo events, taking into account 
the recipe for modelling the momentum spectra of leptons in the rest frame for the semileptonic 
decays of b and c hadrons and for b — > c — £ decays described in a previous paper [|T^]. The 
predicted distribution of Okin is compared to the data in Figure ^ where the contribution of 
events from semileptonic decays of b hadrons is also indicated. Note that the semileptonic 
branching ratios of the individual b hadrons were assumed to be proportional to the lifetimes. 
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Figure 3: The distribution of QeQ2jet for data is shown by the points with error bars. The 
curve indicates the Monte Carlo prediction, assuming Amd = 0.45 ps~^ and Amg = 15 ps~^. 
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Figure 4: The distribution of QeQ2jet for semileptonic decays of neutral B mesons in the Monte 
Carlo. The solid points represent unmixed B mesons while the triangles show the distribution 
for mixed B mesons. 
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5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



«kin 

Figure 5: The distribution of akin for data is shown by the points. The sohd histogram shows 
the Monte Carlo prediction, and the dashed-dotted histogram shows the contribution from 
semileptonic decays of b hadrons. The arrow indicates the selection cut; events to the left of 
the arrow were rejected. 

The overall likelihood of the event sample is taken as 



(15) 



m=l i,j 

where the product is over all the events in the sample. 

6 Results 



In order to obtain information about Am^ and Amg, the likelihood, C, was maximised with 
respect to Am^ for many values of Amg. However, in order to include systematic uncertainties, 
other parameters in the likelihood calculation were also allowed to vary, but under Gaussian 
constraints, i.e. 

^ = -^data X -^constr (1^) 

where £constr is the product of the Gaussian constraints. Thus a multiparameter fit was per- 
formed for each Amg point. This procedure was also employed in reference [^. The parameters 
that were allowed to vary and the constraints on these parameters are listed in Table ^. In 
this table, fs is defined as the production fraction /(b — > B^). The value is constrained both 
by direct measurements, giving a rate of (11.1 ± 2.6)% |15[ relative to all weakly decaying b 
hadrons, and by the measured average mixing rate of b hadrons, x = 0.126 ±0.008 |15| together 



14 



Parameter 


Constraint 


fs 


see text 


fhaxyon 


(13.2 ±4.1)% 




±15% X nominal 


fc 


±20% X nominal 


fuds 


±20% X nominal 




1.549 ± 0.020 ps 




1.03 ±0.06 




1.03 ±0.075 




0.73 ±0.06 


5Qb 


0±0.06 




0±0.06 


6Qb+ 


0±0.06 


'^Qudsc 


0±0.06 




(35± 10)% 



Table 4: Constrained parameters in the fit 



with knowledge of the equivalent parameters, Xd(s), for B°(Bg) mesons and the fraction of b 
baryons, where 

= o f TT^ 1 ' = AmdTd , (17) 

with Td representing the lifetime, and similarly for Xs- This is equivalent to the constraint 
fs = (11.2l|;9)% except that the values of Xd and Xs are calculated from the values of Amj 



and Ams in the fit, together with the appropriate lifetimes, and /baryon, defined as /(b — > b 
baryon), is also taken from the value in the fit. The quoted constraint on /baryon is taken from 
reference |]T3[. When varying /g and /baryon, the following rules are applied: 



/d ± /u ± /s ± /baryon = 1 and fd = /u , (18) 

where fd and /u are the equivalent production fractions for and B"*" respectively. The 
actual fraction of selected leptons coming from B^ mesons (b baryons) is also affected by the 
semileptonic branching ratio, which is assumed to be proportional to the ratio of lifetimes 
Ts/ (Tb) (tav,/ (^b)), where Tg is the Bg lifetime, ta^ the average lifetime for b baryons and (xb) is 
the average lifetime for b hadrons. Similarly, the fraction of selected leptons coming from B"*" 
relative to B[j decays depends on the ratio r^/rj, where r"*" is the B+ lifetime. The parameter 
fh^c^i is a scaling factor affecting the source fractions for b ^ c — > £ processes, sources 6 to 
8. The uncertainty quoted includes uncertainties due to branching fractions, decay modelling 
and detector simulation. The parameter fc is a scale factor multiplying the source fractions 
for sources 2 and 3, and /uds is the corresponding parameter for source 1. The uncertainty 
on fc covers uncertainties due to the charm semileptonic branching ratio, the modelling of the 
akin distribution and lepton background. The uncertainty on f^ds is due to uncertainties in the 
lepton background rates. The lifetime constraints indicated in the table were taken from the 
Particle Data Group review |l^. Note that the existence of two different lifetimes for the Bg 



mass eigenstates is expected to have a negligible effect on this analysis and has been neglected. 
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The parameter SQ^ represents an offset to the Monte Carlo distributions of Q2jet for Z — bb 
events before calculation of rjij. Such an offset moves the distributions for produced b and b 
hadrons closer together or further apart. Thus it causes a change in 77 as a function of Q2jet, and 
is equivalent to an uncertainty in rj of about ±0.02 when 77 = 0.25. Monte Carlo predictions for 
rj were tested to this level in a previous paper |]25[. The parameter (5Qmix has the same definition 



as 6Q\j, but applies only to events where the lepton originates from a mixed B meson. The 
parameter 6Q^+ is also defined in the same way, but applies only to events where the lepton 
originates from B"*" decays. Similarly, 5Qudsc applies only to uu, dd, ss or cc events. The 
parameter /d*. simulates the uncertainty on the rate of b — c — £ decays from B° and B"*" 
mesons relative to all b hadrons due to the effect of varying the rate of decays through a D**. 
A variation in the rate of b — c — ^ £ decays from B^ mesons relative to all b hadrons was also 
considered, but found to have a negligible effect. 

The uncertainty due to the resolution functions, 71 and S, was not taken into account using 
this technique. Instead, these functions were reparametrised without applying the smearing of 
the track impact parameters in the Monte Carlo, and the likelihood analysis repeated using 
the new functions. The effect of fragmentation uncertainties on the resolution functions was 
assessed by using Monte Carlo events generated with {xE)h shifted by 0.02 relative to the central 



value. Such a change represents a shift of over 2a with respect to the measured value but 
is inflated to include the effect of shape uncertainties. Uncertainty due to charm fragmentation 
is expected to have a negligible effect, and was neglected. 



6.1 Determination of Amd 

With Arris fixed to 15 ps^^ the multiparameter fit was performed, with the result 

Amd = 0.444 ± 0.034 ps"^ . 

This error includes systematic components, due to the constraints of Table ^, as well as a 
statistical component. The fitted value of fs was (10.6 ± 1.9)%, and the calculated value of x 
was 0.115. The distribution of t from the fit is superimposed on the data in Figure ^. Shown 
in Figure ^ is the ratio R versus t for |(52jet| > 2, and {QiQ2jet) versus t for |Q2iei| < 2. Here 
R is defined as the fraction of events that have QiQ2jet < 0. The ratio R is not an appropriate 
variable to plot for |(52jet| < 2, because rj changes rapidly with |Q2jet| in this region and, in 
particular, many of the events have |Q2jet| close to 0, where there is no sensitivity. The mean 
value of QeQ2jet is more sensitive than R in this region. Conversely, since t] changes only slowly 
with Q2jet when |Q2jet| > 2, i? is the more appropriate variable in this region. 

The systematic error on Am^ due to tracking resolution was assessed, as described in the 
previous section, to be ±0.004 ps"^. The contribution from the b fragmentation uncertainty 
was found to be ±0.003 ps"^. The effect of the uncertainty in the fraction of B mesons decaying 
to D** on the fitted resolution functions was assessed by repeating the analysis using resolution 
functions determined from a Monte Carlo that did not include D**. The difference between 
the resulting Amd value and the central value quoted above was divided by 3 to correspond 
to an uncertainty of ±10% (a relative error of 30%) in the production of D** per B meson. 
This procedure contributed an uncertainty of ±0.001 ps~^ on Amd. The systematic error due 
to Amg was assessed by taking the largest excursions on the fitted value of Amd when the fit 
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2 4 6 8 10 

t/ps 

Figure 6: The distribution of t for the selected events in the data. The sohd hne represents 
the resuh of the fit, while the dashed line indicates the component from semileptonic decays of 
b hadrons, sources 9, 10 and 11. The insert shows the same distribution on a linear scale. 



17 




Figure 7: a) The ratio R versus t for data with |(52jet| > 2. b) (.QnQijet) versus t for data with 
\Q2jei\ < 2. The hne shows the result of the fit when Arris = 15 ps^^ in both a) and b). 



was repeated for different values of Arris in the range Arris > 6.6 ps ^. We consider the range 
Arris < 6.6 ps^^ to be excluded by previous measurements |2^, |2^, |2^, |2^. The error on Am^ 



from this source was lo;ooo2 ps~^. The errors on Arrid resulting from different ranges of Arrig 
are given in Table ||. The uncertainty due to the weighting function w{\Q2jet\) was found to be 
negligible. 



Lower limit on Arris (ps ^) 


Systematic uncertainty on Amd(ps ^) 


3 


+0.014 -0.0002 


6.6 


+0.009 -0.0002 


8 


+0.004 -0.0002 


10 


+0.001 -0.0002 



Table 5: The systematic error on Am^ resulting from different ranges of Arris values. 



The final result for Am^ is Am^ = 0.444 + 0.034 ps~^, where, as above, the first error 
has a systematic component. To aid comparison with previous measurements, the systematic 
components to the fit error (due to the constraints of Table ^ were calculated by determining 
the effect of changing each constraint by the error given in Table one at a time. The results 
are shown in Table We can now rewrite the result for Am^ as 

Amd = 0.444 + 0.029 tof^ Ps"^ , 

where the first error is essentially statistical and the second is systematic. Note that all con- 
strained parameters were consistent with the constraints to within one standard deviation, 
except /b^c^^, /uds and SQudsc where the fitted values were larger than the constraint values, 
but were consistent at the level of two standard deviations. 
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Origin 


Systematic error on Awd in ps ^ 


fs 

J baryon 


-U.Uii 
+0.010 
+0.005 
-0.004 


Jb->c-»^ 




f 
/c 


4-n nnn 


J nas 


±0.000 




±0.000 




±0.004 


Ts/Td 


^0.002 

+0.007 
-0.006 




^0.001 


'^Qmix 


^0.007 




±0.006 




^0.003 


/d** 


±0.000 


Tracking resolution 


±0.004 


b fragmentation 


±0.003 


Resolution effect of D** 


^0.001 


Ams variation 


+0.009 
-0.000 


Total 


+U.U2U 
-0.017 



Table 6: Summary of systematic errors on Amd 
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6.2 Study of ArUs 



In order to constrain Amg, we consider two different techniques. For the first technique, we plot 
in Figure |^ the value of —A log£ as a function of Auis, where A log£ indicates the difference 
in log£ relative to the maximum value of log£. The solid points include the effects of the 
systematic uncertainties. They were obtained by allowing all parameters, including Amd, to 
vary at each value of Arris. To account for systematic errors due to the description of the 
resolution, four different sets of resolution functions were used. In addition to the default 
functions, we used the resolution functions obtained without the impact parameter smearing 
in the Monte Carlo, and those obtained using the Monte Carlo with different fragmentation 
parameters. A fourth choice of resolution function was obtained by using only semileptonic 
Bg decays to parametrise sources 5,9,10 and 11 (using the nominal tracking resolution). For 
each of the four sets of resolution functions the full minimisation of the other parameters was 
performed at each Amg point, giving four distributions of Alog£, with respect to the point of 
maximum C relevant to each set, versus Arris. The solid points were obtained simply by taking 
the minimum value of Alog£ from the four choices at each Arrig point. As a result of this 
procedure, the points for Arris < 1 ps^^ and for Arris > 3.6 ps^^ are taken from the resolution 
functions without the impact parameter smearing, while those for 1 ps^^ < Arrig < 3.6 ps~^ 
come from the resolution functions obtained with the modified fragmentation parameters. The 
value of Alog£ at Arris = was found to be 2.6 (including systematic errors). A limit on 
Arrig is extracted by finding the intersection of these points with a limit curve calculated using 
a Monte Carlo technique. This technique consisted of performing fast simulations of the data 
sample, assuming various values of Amg. The simulated samples were then analysed in the 
same way as the data (including the multiparameter fit), and a distribution of Alog£ obtained 
at each Amg point. The limit curve is defined by the value of Alog£ above which lie only 5% 
of the simulated samples. A total of 6000 samples were simulated, distributed over four values 
of Arris. Three of the limit points and the limit curve, obtained by linear interpolation, are 
indicated in Figure The fourth limit point is at Arris = 50 ps~^. 

We extract a limit Arris > 3.1 ps^^ at 95% confidence level. The region 5.0 ps^^ < Arris < 
7.6 ps~^ is also excluded. The dotted curve shows the result of performing the full minimisation 
at the global minimum. Arris = 3.6 ps~^, using the default resolution functions, and then 
scanning through Arris while keeping all other parameters fixed. The fitted value of fs, of 
particular relevance for ArUs measurements, was (11.5 ± 1.9)%. 



A second technique for studying Arris, known as the amplitude method was also used. 
The method consists of replacing the cos Arrist' term in the physics function with A cos Arrist', 
where A is known as the amplitude. The parameter A may be fitted at each value of Amg, and 
if the value is smaller than unity then the value of Arris may be excluded with a confidence that 
depends only on the value and uncertainty of A. The advantage of such a technique over an 
exclusion based on A log C is that A is measured with errors that are approximately Gaussian. 
This makes it easy to combine results and compare sensitivities from different analyses or 
experiments. The calculation of a limit is also relatively straightforward. 

For this analysis, A was simply added as an extra parameter to the fit. All systematic 
uncertainties were handled in the same way as for the Ama result. The fitted value of A as a 
function of Arris is shown in Figure To facilitate combination with other results, the central 
values and a breakdown of the error contributions in steps of 1 ps~^ are documented in the 
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Anig / ps' 

Figure 8: The dotted line shows the distribution of — Alog£ with respect to the maximum 
value of logC as a function of ArUs, while keeping all other parameters fixed to values optimised 
at ArUs — 3.63 ps~^. The solid points include the effects of systematic uncertainties. The open 
triangles represent the position of the 95% confidence level, and the horizontal hne joins together 
these points, increased by one standard deviation. The hatched region indicates values of Amg 
excluded at 95% confidence level by this analysis. 
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Appendix. The results are qualitatively in agreement with the A log L method. There is a 
possible signal near Amg = 4 ps~^, as seen also in the Alog£ plot, where A is close to 1 and 
significantly separated from 0. Previous measurements ^ exclude this region with a high 
confidence, so that this excursion is more likely to be a statistical fluctuation. To determine 
exclusion regions at 95% confidence level, at a given value of Amg we represent the measured 
value of y4 as a Gaussian distribution function G{A — \i, a a), where [i is the central value and oa 
is the measurement error. Two alternative methods are then considered to determine whether 
the value of Amo is excluded: 



a) values are excluded where the probability of measuring an amplitude lower than that ob- 

served is less than 5%, if that value of Amg was the correct one, i.e. 

/oo 
G(A-^,aA)dA < 0.05 , (19) 

or 

b) the same definition, but limited to the positive region, i.e. 



/o~G(A-/i,a^)dA 



< 0.05 . (20) 



The first definition gives a true 95% confidence level, in the sense that there is a 5% probability 
to exclude the true value. However, it is not protected against setting hmits well beyond the 
experimental sensitivity. The second definition makes use of the fact that the predicted value 
of A lies between and 1, regardless of the value of Arris . It is automatically protected against 
setting limits beyond the sensitivity. For a true value of Amg well beyond the sensitivity, 
method a) would exclude the true value in 5% of the experiments, while for method b) this 
percentage would tend towards zero. 

For method a), excluded regions are defined by A + 1.645cryi < 1, where a a represents 
the total error on A. This gives the result Amg > 2.9 ps~^ at 95% confidence level. The 
small region 6.4 ps~^ < Amg < 6.7 ps~^ is also excluded. For method b), we find a limit 
Amg > 2.9 ps~^, and higher Amg regions are not excluded. The difference between the two 
lower limits is 0.03 ps^^. Both results are similar to those obtained from the A log£ technique. 

We define a measure of the sensitivity of this experiment as the value of Amg that would 
be excluded by method a) if A were measured to be 0, i.e. where 1.645(7^1 = 1. The sensitivity 
is 4.8 ps^^ 



6.3 Data tests 

The plots of R (the fraction of events with QiQ^jet < 0) versus t and {QeQ2jet) versus t shown 
in Figure |^ are repeated in Figure |10|, but restricted to the time windows —0.25 ps < t < 3 ps 
and —0.5 ps < t < 3 ps respectively. Also shown in these plots are the fitted curves assuming 
Amg = 3.6 ps~^ and Amg = 15 ps~^ using the nominal tracking resolution. Calculating the 
summed over both plots (over the full time range) yields values of 44.8 for Amg = 3.6 ps~^ 
and 48.0 for Amg = 15 ps~^ for 39 bins in each case. The difference in log£ between these 
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Figure 9: The result of the amphtude fit. The points show the fitted value of the amplitude, 
^4, as a function of Amg, and the error bars include all systematic errors. 
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Figure 10: a) The ratio R versus t for data with \Q2jet\ > 2, restricted to the time window 
—0.25 ps < t < 3 ps. b) The value of {QeQ2jet) versus t for data with |Q2iet| < 2, for the time 
window —0.5 ps < t < 3 ps. The sohd hne shows the result of the fit when Amg = 3.6 ps~^, 
and the dashed-dotted hne corresponds to Arus = 15 ps~^. 

two values of Amg is 2.5 (again using the nominal tracking resolution), showing a behaviour 
similar to that of the calculated values. This test suggests that the global minimum at 
Ams = 3.6 ps~^ is not an artefact of the fitting procedure, but is rather a property of the data. 

To test the stability of the results, the Ama fit was repeated using a tighter selection: the 
ctkin cut was changed from 0.7 to 0.85. The fit result was 0.438 ± 0.039 ps"\ consistent with 
the previous result of 0.444 ± 0.034 ps~^ 

6.4 Monte Carlo tests 

The multiparameter constrained fit described above was performed on a Monte Carlo sample 
based on four million hadronic events generated with Amj = 0.438 ps~^ and an infinite value of 
Ams, with Ams fixed to a large value in the fit. The fitted value for Amj was 0.430±0.030 ps~^. 
The fitted values for three of the parameters, (rb), /uds and /b^c^^j were more than 2a from 
their true values. If the uncertainty due to the resolution functions and their parametrisation 
is taken into account, then all values of all parameters are consistent with their true values. A 
possible bias associated with these parameters was investigated in the data by repeating the fit 
with all three of these parameters fixed to their nominal values. The result for Ama was shifted 
by +0.003 ps~^. The effect on the fitted amplitude for oscillations at Arris = 3.75 ps~^ was a 
shift of 0.18, which is easily covered by the total error of 0.50. The possible bias does not cause 
a significant problem for the results, and is addressed by resolution function uncertainties. 
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To test the sensitivity to Amg, the single sample of four million simulated events was used 
to simulate event samples with Am^ = 0.45 ps~^ and a range of values for Amg: 1 ps~^, 2 ps~^, 
3 ps~^, 4 ps~^, 5 ps~^, 6 ps~^, 8 ps~^, 10 ps~^ and 15 ps~^ The values of the fit parameters 
other than Amg were fixed at their known values in the simulation, and log£ calculated as a 
function of Arris for each simulated event sample. The resulting plots of — Alog£ versus Arris 



(relative to the point of maximum likelihood) for simulated data are shown in Figure 11. It can 
be seen that A log£, calculated at the generated value, is in most cases smaller than 2.8, which 
corresponds roughly to a 95% confidence level. The exceptions are the samples generated with 
Arris = 8 ps^\ where the generated value is just excluded, and Arris = 15 ps~\ but note that 
the samples are statistically correlated with each other. It is interesting to note that in the 
case of the sample generated with Arris = 15 ps~^, an apparent signal would have been seen at 
a low value of Arris. Such cases are properly included in the fast Monte Carlo simulations. 

Using the same simulated samples, the equivalent test was performed for the amplitude 
method. All fit parameters were fixed with the exception of A, which was fitted as a function 
of Amg, for each sample. The results are shown in Figure |12[ For the sample generated with 
Amg = 8 ps~^, the generated value of Amg would be excluded by method a), but not by method 
b). In none of the other samples would the generated Amg value be excluded by either method. 



7 Test of CP(T) conservation 

The multiparameter fit was repeated with two extra free parameters, Reee and Im^B, with 
the physics function modified as described in section 5. The results were first tested on four 
Monte Carlo samples simulated with and without CP and CPT violation. The four Monte Carlo 
samples are statistically correlated and are each equivalent to four million hadronic events. The 
results of this test, in which Amg was fixed to a large value (and A fixed to 1), but all other 
parameters were free to vary, are given in Table |^. 



Generated Re ee 


Generated Im^B 


Fitted RecB 


Fitted ImSB 





0.1 
0.2 




0.1 




-0.014 ±0.009 
-0.015 ±0.010 
0.077 tril 

n 1 77 +0.018 
U.J- ' ' -0.016 


0.005 ±0.014 
0.096 ±0.015 
-0.002 iril 
-0.013 ±0.016 



Table 7: Results of fits for Reee and Im^B using simulated data. There is a strong statistical 
correlation between the different samples. 



The result of the fit to the data, with Amg fixed to 15 ps , was 



Re ee 
Im 6b 



-0.006 ±0.010 
-0.020 ±0.017 



The distribution of -Riept, the fraction of leptons that are negatively charged, versus recon- 



structed time is shown in Figure |T3| for the data, with the fit result superimposed. In this 
figure, the data is divided into four categories according to the sign of QeQ2jet and whether 
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Figure 11: The value of — Alog£, relative to the maximum value of log£. as a function of 
Arris for data simulated using nine different Arris values, indicated by the dashed line in each 
case. All other parameters in the hkelihood calculation were fixed. The Monte Carlo event 
statistics are equivalent to 4 million hadronic Z° decays. Note that the different values of Arris 
are not statistically independent. Horizontal lines are drawn on each plot at — Alog>C = 2.8, 
corresponding roughly to a 95% confidence level. 
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Figure 12: The fitted value of A as a function of ArUs for data simulated using nine different 
Arris values, indicated by the arrow in each case. All parameters other than A were fixed in 
the fit. The Monte Carlo event statistics are equivalent to 4 million hadronic Z° decays, but 
the nine samples are statistically correlated. 
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|Q2jei| is smaller or larger than 2. Note that, following equation 10, the data with QiQ2jet < 
(enhanced in mixed events) are more sensitive to Reee, while the data with QiQ2jet > are 
more sensitive to Im^B- For comparison, the predicted curves assuming Reee = 0.1, Im^B = 
and also Reee = 0, Im^B = 0.1 are included in the figure. These curves were obtained keeping 
all other parameters fixed at their fitted values. 

As in the case of the Am^ analysis, the errors quoted contain both statistical and system- 
atic components. The systematic errors are shown in Table |^, including the errors from the 
resolution function, calculated as in the Arud analysis. The uncertainty due to the weighting 
function w{\Q2jet\) is also indicated, estimated by shifting the jet charge by 0.06 for Monte Carlo 
simulated with x = 0.5 in the opposite jet. Such a shift causes a 10% change in w{\Q2jet\)- 





11/11 ui uli rvc eg 


11/1 lUi Uli llli eg 


1-/1 1 Ul Uli IXctg, 1111 eg — u 


/s 


±0.000 


±0.000 


±0.000 


/baryon 


±0.002 


±0.002 


±0.000 




±0.000 


±0.000 


±0.000 


fc 


±0.000 


±0.000 


±0.000 


fuds 


±0.000 


^0.001 


±0.000 


(rb) 


TO.001 


^0.001 


±0.000 


r+/rd 


±0.000 


±0.001 


±0.000 


rs/rd 


^0.001 


±0.000 


±0.000 




^0.001 


^0.001 


±0.000 


(^Qb 


±0.000 


±0.000 


±0.000 




±0.001 


±0.001 


±0.000 


6Qb+ 


±0.001 


±0.002 


±0.000 




±0.001 


±0.001 


±0.000 


/d** 


±0.000 


±0.000 


±0.000 


Tracking resolution 


±0.001 


±0.001 


±0.000 


b fragmentation 


±0.001 


±0.002 


±0.000 


Resolution effect of D** 


±0.000 


±0.000 


±0.000 


Arris variation 


±0.000 


±0.000 


±0.000 


w{\Q2jet\) 


±0.001 


±0.000 


±0.000 


Jet charge asymmetry 


±0.004 


±0.003 


±0.003 


Lepton charge asymmetry 


±0.001 


±0.001 


±0.001 


Momentum asymmetry 


±0.001 


±0.001 


±0.001 


Lepton background 


±0.002 


±0.003 


±0.001 


cos 6 asymmetry 


±0.000 


±0.000 


±0.000 


Total 


±0.006 


±0.006 


±0.003 



Table 8: Summary of systematic errors on Reee and Im^B- The fourth column gives the 
systematic errors on Reee when Im^B is fixed to 0. 



In addition to these uncertainties, extra sources of error that affect only the CP(T) param- 
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Figure 13: The distribution of i?iept versus reconstructed time for the data, with the result 

of the CPT fit superimposed. The data is divided into four categories : a) QtQ2jet > and 
\Q2jet\ > 2, b) QiQ2jet < and \Q2jet\ > 2, c) QtQ2jet > and \Q2jet\ < 2, d) QiQ2jet < and 
|(52jet| < 2. Also shown are the predicted curves for Reee = 0.1, Im^B = and for Reee = 0, 
Im^B = 0.1. 
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eters had to be considered. Firstly, the mean jet charge for produced b quarks is not exactly 
opposite that for produced b quarks. This is caused by a slight asymmetry in the jet chamber 
in the efficiency to detect low momentum positively and negatively charged tracks. The asym- 
metry arises from the geometry of the detector and is therefore simulated in the Monte Carlo. 
The mean jet charge, averaged over b and b quarks, is found to be +0.020 ± 0.006, where the 
error is statistical only. This result was checked by studying the charge asymmetry for tracks 
in hadronic events for both the data and the Monte Carlo. The mean jet charge could then 
be calculated by weighting the charge asymmetry appropriately as a function of momentum, 
according to the contribution to the jet charge. Consistent results were obtained whether the 
weighting was determined from data or Monte Carlo. The resulting mean jet charge values were 
0.008 for data and 0.013 for Monte Carlo. This calculated Monte Carlo result is consistent with 
the value directly observed in the Monte Carlo. The final value for the data is determined by 
scaling 0.008 by 0.020/0.013, assigning an error of 0.006 due to Monte Carlo statistics, and 
0.008 due to the discrepancy seen between Monte Carlo and data, to give 0.012 ± 0.010. This 
tracking asymmetry could also have a minor effect on the efficiency for selecting leptons (re- 
quired to have p > 2 GeV/c). In this case, the statistics of the Monte Carlo are insufficient to 
investigate the effect directly. The effect was checked by studying inclusive tracks, passing the 
same track quality cuts as the leptons, as a function of p. A possible asymmetry of 6 x 10~^ 
was deduced for a value of (l/p) = 0.14 GeV~^, corresponding to the lepton sample. This was 
assigned as a systematic uncertainty. Another cause of uncertainty would be an asymmetry 
in the momentum measurements for positive and negative tracks, which could be caused by a 
slight rotation of the outside of jet chamber relative to the inside. This effect was investigated 
by studying Z° — events in the data. A difference in (l/p) of (7± 2) x 10~^ GeV~^ was 
observed between /i~ and /i"*" particles. Such a shift could cause a lepton selection asymmetry 
of 5 X 10"'*, assigned as a systematic error. 

The results are also sensitive to a charge asymmetry in the lepton background. This asym- 
metry was studied in the data using electron and muon candidates that just failed the lepton 
selection criteria. The charge asymmetry was found to be (1.9 ± 0.7)%, with the largest com- 
ponent due to muon candidates resulting from non-interacting kaons and kaon punchthrough. 
The central values quoted above were corrected for this effect. 

Finally, a cos 6' asymmetry for lepton identification in the detector would cause a charge 
asymmetry through the Z ^ bb forward-backward asymmetry, measured by OPAL to be (9.06± 
0.51 ± 0.23)% at the peak |3^. CP(T) violation would not induce a cos 6' asymmetry. A 



lepton cos6' asymmetry of (— 1.7± 0.3)% (more leptons at cos 6* < 0) was observed in the data. 
This can be understood as due to an observed average shift of about +0.5 cm of the beam spot 
along the z-axis relative to the centre of the detector. The central values of Re ee and Im 6b 
quoted above were adjusted by —0.002 and —0.001, respectively, for this asymmetry and 20% 
(representing the relative error on the cos 6 asymmetry) of the shifts taken as systematic error. 

Splitting the fit errors into their statistical and systematic components and including the 
extra systematic errors described above, the final results are : 

ReeB = -0.006 + 0.010 + 0.006 
Im Sb = -0.020 + 0.016 + 0.006 . 

These results neglect CP violation in the system. The sensitivity of the results to Reeeo, 
where eeo is e for the system, was gauged by repeating the fit assuming a large CP asymmetry. 
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ReeBO = —0.05. The value of Reee was shifted by +0.015, while Iiii^b was shifted by +0.006. 
Note that the influence of Im 6-^0 7^ on the results is negligible, because the effect is smeared 
out by the rapid oscillations (see equation 10). 

CPT is usually assumed to be a good symmetry, i.e. = 0. In the results presented above, 
there is a correlation of order +30% between the Reee and Im^B results, so we can obtain 
more precise Reee results if Im^B is set to 0. The fit was repeated with Im^B = to give 

Reee = 0.002 + 0.007 , 

where the fit contains both statistical and systematic errors. The systematic error calculation 
was performed as described above, and the errors are included in Table ^. The final result is 

Re eg = 0.002 + 0.007 + 0.003 . 

In this case, the effect of Reeeo = —0.05 causes a shift of +0.013 in Reee. 



8 Conclusions 

A sample of inclusive lepton events was used to study B[| and oscillations. An estimate of 
the decay proper time of the inferred b hadron was reconstructed for each lepton candidate, 
and a jet charge technique was used to tag the produced B flavour. We measure 

Amd = 0.444 + 0.029 tofij ps"^ . 



This result is consistent with previous measurements [24, 25, pTl, and represents the most 



precise result for a single technique. Taking account of the common systematic errors, we 



combine this measurement with previous OPAL results |24, 25 1 to give 



n A7K+0-022 +0.018 



The small statistical correlations between the results have a negligible effect. 

By studying the behaviour of log £ as a function of Amg, we are able to exclude the regions 
Arris < 3.1 ps~^ and 5.0 ps~^ < Arrig < 7.6 ps~^ at 95% confidence level. Using an amplitude 
method instead, the lower limit of 3.1 ps^^ is slightly weakened to 2.9 ps~^, and the region 
5.0 ps^^ < Arris < 7.6 ps~^ is either only partially excluded or not excluded at all, depending 
on the details of the confidence level calculation. The lower limit that would be obtained using 
the amplitude method, were the amplitude measured to be zero, is 4.8 ps~^, a measure of the 
experimental sensitivity. This result is consistent with previous results p4|, |26|, |2^, ^ , of which 
the most constraining 



29| quotes a lower limit of 6.6 ps ^ at 95% confidence level. 



By studying the charge symmetry of the B[} mixing structure, we are able to constrain 
possible CP and CPT violating effects. We measure simultaneously the indirect CP violation 
parameter 

ReeB = -0.006 + 0.010 + 0.006 
and the indirect CPT violation parameter 



Im^B = -0.020 + 0.016 + 0.006 
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If alternatively we invoke CPT symmetry, then we obtain 



ReeB = 0.002 ± 0.007 ± 0.003 



The ReeB measurement, whether or not CPT symmetry is invoked, is consistent with, but 
more precise than, previous measurements from CLEO |^ and CDF [Q. The Im^B result 
represents the first measurement of this quantity in the B system. 



Appendix: Amplitude results 

Section 6.2 describes the amplitude results for the Arris study. We detail in Table ^ the central 
values of the amplitude and the breakdown of the error contributions in steps of 1 ps~^. This 
information is essential for a correct combination of the results from this paper with other anal- 
yses. The systematic uncertainties for intermediate points may be estimated by interpolation. 



Acknowledgements 



We particularly wish to thank the SL Division for the efficient operation of the LEP acceler- 
ator and for their continuing close cooperation with our experimental group. We thank our 
colleagues from CEA, DAPNIA/SPP, CE-Saclay for their efforts over the years on the time- 
of-fiight and trigger systems which we continue to use. In addition to the support staff at our 
own institutions we are pleased to acknowledge the 
Department of Energy, USA, 
National Science Foundation, USA, 
Particle Physics and Astronomy Research Council, UK, 
Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 

Minerva Gesellschaft, 

Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Science and Culture (the Monbusho) and a grant under the 
Monbusho International Science Research Program, 
German Israeli Bi-national Science Foundation (GIF), 

Direction des Sciences de la Matiere du Commissariat a I'Energie Atomique, France, 
Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologic, Germany, 
National Research Council of Canada, 

Hungarian Foundation for Scientific Research, OTKA T-016660, T023793 and OTKA F-023259. 



32 



Arris (ps ^) 





1 


2 


3 


4 


5 


6 


7 


A 


0.53 


0.29 


-0.32 


0.44 


0.99 


0.60 


0.16 


-0.35 


_.stat 


±0.18 


±0.19 


±0.29 


±0.35 


±0.46 


±0.59 


±0.75 


±0.91 


fs 


±0.09 


±0.11 


±0.01 


^0.08 


TO.23 


to.io 


TO.08 


±0.00 


/baryon 


±0.01 


±0.01 


TO.04 


±0.02 


=F0.02 


±0.00 


±0.02 


±0.03 




±0.04 


±0.00 


TO.02 


TO.02 


TO.01 


TO.02 


TO.01 


±0.00 




±0.01 


±0.00 


±0.00 


±0.01 


±0.04 


±0.01 


±0.02 


±0.03 


fc 


±0.00 


^0.01 


^0.01 


±0.00 


^0.01 


^0.01 


TO.01 


±0.00 


/uds 


±0.00 


±0.01 


±0.01 


±0.00 


±0.01 


±0.01 


±0.01 


±0.00 


(rb) 


±0.00 


±0.01 


±0.00 


±0.00 


TO.01 


±0.00 


±0.00 


TO.01 


r+/rd 


::F0.03 


TO.03 


±0.01 


to.oi 


to.oi 


TO.OI 


TO.02 


±0.00 




±0.02 


±0.02 


±0.00 


TO.03 


TO.01 


TO.02 


TO.01 


TO.03 




^0.01 


±0.01 


^0.02 


±0.00 


^0.02 


±0.00 


±0.00 


±0.01 


•^Qb 


TO.06 


TO.02 


±0.02 


TO.02 


TO.06 


±0.00 


TO.01 


±0.05 




TO.01 


TO.01 


±0.00 


TO.02 


TO.08 


TO.03 


TO.02 


±0.01 




TO.03 


TO.02 


TO.02 


±0.02 


±0.04 


±0.01 


±0.01 


TO.04 


<^Qudsc 


^0.01 


±0.00 


^0.05 


:F0.07 


q:0.09 


to.io 


TO.10 


T0.11 




±0.01 


±0.01 


±0.00 


±0.00 


±0.04 


±0.01 


TO.01 


TO.01 


Tracking resolution 


±0.01 


±0.03 


±0.08 


±0.11 


±0.15 


±0.14 


±0.12 


±0.01 


b fragmentation 


±0.00 


±0.02 


±0.02 


±0.02 


±0.01 


±0.00 


±0.01 


±0.04 


Bg resolution function 


±0.01 


±0.03 


±0.03 


±0.01 


±0.04 


±0.09 


±0.01 


±0.06 


syst 


±0.13 


±0.14 


±0.12 


±0.16 


±0.31 


±0.23 


±0.18 


±0.16 



Arris (ps ^) 


8 


9 


10 


11 


12 


13 


14 


15 


A 


1.11 


1.27 


3.14 


2.32 


2.36 


3.45 


5.23 


8.90 


'^A 


±1.06 


±1.25 


±1.43 


±1.69 


±1.97 


±2.26 


±2.64 


±3.01 


fs 


TO.18 


TO.04 


TO.46 


TO.09 


TO.49 


±0.00 


TO.71 


Tl.26 


/baryon 


TO.03 


TO.01 


TO.12 


±0.03 


TO.07 


±0.00 


TO.09 


TO.41 




±0.00 


TO.03 


TO.04 


TO.01 


TO.01 


TO.04 


TO.12 


TO.14 




±0.03 


±0.02 


±0.03 


±0.01 


±0.01 


±0.01 


±0.01 


±0.05 


fc 


TO.03 


TO.02 


TO.07 


TO.04 


TO.04 


TO.03 


TO.13 


TO.23 


/uds 


±0.03 


±0.01 


±0.07 


±0.03 


±0.04 


±0.03 


±0.13 


±0.22 




TO.01 


TO.01 


±0.01 


±0.00 


±0.02 


TO.01 


±0.00 


±0.01 


r+/rd 


±0.00 


TO.04 


±0.00 


±0.01 


TO.02 


TO.01 


±0.02 


±0.02 


Ts/Td 


±0.00 


TO.09 


±0.04 


TO.03 


TO.01 


TO.04 


to.io 


±0.02 


'^Ab/rd 


TO.01 


TO.02 


±0.00 


±0.02 


±0.04 


±0.00 


±0.04 


TO.07 




TO.01 


±0.01 


TO.07 


±0.02 


TO.10 


±0.02 


T0.11 


TO.33 




TO.06 


TO.03 


TO.12 


TO.03 


TO.13 


TO.03 


TO.18 


TO.38 




±0.01 


±0.02 


±0.05 


±0.04 


±0.08 


±0.02 


±0.03 


±0.23 




TO.12 


TO.14 


TO.14 


TO.17 


TO.20 


TO.19 


TO.21 


TO.23 


/d" 


±0.01 


±0.01 


±0.00 


TO.03 


TO.02 


±0.02 


±0.04 


±0.04 


Tracking resolution 


±0.02 


±0.01 


±0.06 


±0.06 


±0.18 


±0.29 


±0.45 


±0.61 


b fragmentation 


±0.38 


±0.47 


±0.31 


±0.10 


±0.45 


±0.79 


±0.61 


±0.43 


Bg resolution function 


±0.31 


±0.45 


±0.38 


±0.52 


±0.55 


±0.59 


±1.07 


±1.56 


syst 


±0.55 


±0.68 


±0.72 


±0.58 


±0.93 


±1.05 


±1.54 


±2.29 



Table 9: Amplitude results with the breakdown of systematic errors. 
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